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ABSTRACT 

We study the structure of an outer-magnetospheric gap around a rotating neutron 
star. Migratory electrons and positrons are accelerated by the electric field exerted in the 
gap and radiate copious gamma-rays via curvature process. Some of the gamma-rays 
materialize by colliding with the X-rays illuminating the gap, leading to a pair pro- 
duction cascade. The replenished charges partially screen the original acceleration field, 
which is self-consistently solved from the Poisson equation, together with the Boltzmann 
equations for gamma-rays and the continuity equations for particles. We demonstrate 
that it is difficult to detect the TeV emission due to Compton upscatterings in the gap, 
by the current ground-based telescopes. 

Subject headings: gamma-rays: observation - gamma-rays: theory - magnetic field - 
X-rays: observation 
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1. Introduction 



2. Basic Equations and Boundary Conditions 



In the very high energy (VHE) region above 100 
GeV, positive detections of radiation at a high confi- 
dence level have been reported from the direction of 
the Crab, B1706-44, and Vela pulsars (Bowden et al. 
1993; Nel et al. 1993; Edwards et al. 1994; Yoshikoshi 
et al. 1997; see also Kifune 1996 for a review), by 
virtue of the technique of imaging Cerenkov light from 
extensive air showers. However, as for pulsed TcV 
radiation, only the upper limits have been obtained 
from these pulsars. If the VHE emission originates the 
pulsar magnetosphere, rather than the extended neb- 
ula, significant fractions of them are expected to show 
a pulsation. Therefore, the lack of pulsed TeV emis- 
sions provides a severe constraint on the modeling of 
particle acceleration zones in pulsar magnetospheres. 

In the picture of Cheng, Ho, and Ruderman (1986a, b), 
who developed a version of an outer magnetosphere 
7-ray emission zone, the magnetosphere should be op- 
tically thick for pair production in order that the TeV 
flux is absorbed to be unobscrvable. This, in turn, 
requires very high luminosities of infrared photons, 
which are generally orders of magnitude larger than 
the observed values (Usov 1994). We are therefore 
motivated by the need to contrive an outer gap model 
which produces less TeV emission with a moderate in- 
frared luminosity. 

More recently Romani (1996) showed that the 
pulsed TeV flux can be somewhat less than 1% of 
the pulsed GeV flux, assuming the functional form of 
the acceleration field as E\\ cx 1/r, where r denotes 
the distance from the star. Subsequently, Hirotani 
(1999, Paper IV) revealed semi-analytically that the 
TeV fluxes from outer gaps are unobservable, by con- 
sidering a 'gap closure condition' so that a single pair 
produces copious 7-rays one of which materialize as 
a pair on average (see also Zhang & Cheng 1998 for 
another kind of closure condition). In this letter, in- 
stead of imposing the closure condition, we investigate 
the spatial distribution of £j|, by solving the Pois- 
son equation and the Boltzmann equations simulta- 
neously. This method was first developed by Beskin 
et al. (1992) for black-hole magnetospheres and ap- 
plied to pulsar magnetospheres by Hirotani & Shibata 
(1999a,b,c; hereafter Paper I, II, III). 

We formulate basic equations in § 2, demonstrate 
the invisibility of pulsed TeV fluxes in § 3, and discuss 
the validity of assumptions in § 4. 



2.1. Reduction of Poisson Equation 

To simplify the geometry, let us introduce a recti- 
linear coordinates for a region around the null surface. 
Suppose that the magnetic field lines are straight lines 
parallel to the x axis, where x increases outwardly. 
We approximate that the null surface (x = 0) is per- 
pendicular to the x axis. It is convenient to non- 
dimcnsionalize the length scales by a typical Debey 
scale length c/oj p (see the discussion in §3.1 in Pa- 
per I for details), where 
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e refers to the magnitude of the charge on the electron, 
SI2 the angular frequency (SI) of the neutron star in 
10 2 rad s _1 unit, c the velocity of light, and B$ the 
magnetic field strength (B) at x = in 10 5 G unit. 
Then the dimcnsionless coordinate variable becomes 

£ = ^ x = 6.25 x 10- 4 fi 2 1/2 B5 1/2 *. (2) 
c 

Assuming the trans-field thickness is comparable 
to or greater than the longitudinal width, we Fourier- 
analyze the Poisson equation in the trans-field direc- 
tions to obtain (Paper I, II, III) 



(3) 



■n+(0-n_(0 



,2 V 



d£ A A 

where A± = (uj p /c)D± and <p(£) = e$(x)/(m c c' 
D±_ is a typical gap thickness in the trans-field direc- 
tions and $ the electrostatic potential. The particle 
densities are normalized in terms of a typical value of 
the Goldrcich- Julian density as 
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where N±(x) are the spatial number density of e ±5 s. 
We evaluate it with the Newtonian value as 
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where a; refers to the inclination angle of the mag- 
netic moment, 9 the clatitude angle of the position at 
which pgj is measured, w the distance from the ro- 
tation axis, and wlc = c/Sl the light-cylinder radius. 



2 



2.2. Boltzmann equations 

Let us first introduce the following dimensionlcss 
7-ray densities in the dimensionless energy interval 
between fa_i and fa: 
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where G±(x, e 7 ) are the distribution functions of 7- 
ray photons having momentum ±rn c ce 7 In this letter, 
we set fa — 10; therefore, the lowest 7-ray energy 
considered is 10m c c 2 . In this letter, we divide the 7- 
ray spectra into 13 energy bins as follows: fa — 30, 
fax = 100, fa = 200, fa = 300, fa = 600, fa = 10 3 , 
ft = 2x 10 3 , fa = 3 x 10 3 , ft=6x 10 3 , fa a = 10 4 , 
fti = 2x 10 4 , /3 12 = 3 x 10 4 , /? 13 = 10 5 . 

Assuming for simplicity that both electrostatic and 
curvature-radiation-reaction forces cancel out each 
other, we obtain the following continuity equations 
for particles (eqs. [31] and [32] in Paper III): 
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The factor 1 — q denotes the projection effect due to 
the three-dimensional motion of the particles, which 
was not considered in Papers I, II, and III. Since the 
electric field is assumed to be positive in the gap, e +, s 
(or e~'s) migrate outwards (or inwards). 

The X-ray field affects the gap structure through 
the pair-production rate r]p±. Putting ej = (fa + 
fa_i)/2, we can evaluate rjp± as 
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where dN x /de x denotes the number density of the 
X-rays between energies e x and e x + de x ; o~ p is the 
pair production cross section when two photons col- 
lide with an angle cos -1 fi c (Berestetskii et al. 1989; 
eqs. [21]-[23] in Paper II). For surface blackbody X- 
rays, we obtain ^ c = cos a bb sin C , where </> a bb = 
tan _1 (ci7/ct7Lc) is the aberration angle. For power- 
law, magnetospheric X-rays, we assume fj, c — cos(0 .5W / zulc] 
where W is the full gap width, for both inwardly 
and outwardly propagating 7-rays. As the lower 
bound of the integral, the threshold energy, e-th = 



A combination of the two continuity equations (g|) 
gives the current density jo = n+(£) + which 
is constant for £. When j = 1.0, the current density 
equals the typical Goldreich- Julian value, HB/(2tt). 

We next consider the Boltzmann equations for 7- 
rays. Integrating the 7-ray Boltzmann equations be- 
tween fa-i and fa, we obtain (eq. [39] in Paper III) 
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where i = 1,2, • • •,m(= 13). The 7-ray production 
rate via curvature radiation, 77c 1 , are defined by equa- 
tion (40) in Paper III. The particle motion is assumed 
to be mono-energetic and the Lorentz factor saturates 
in a balance between electrostatic acceleration and 
the curvature drag force (see eq. [19] in Paper III). 

2.3. Boundary Conditions 

In the same manner as in Papers I, II, and III, we 
impose the boundary conditions at the inner (£ = £1) 
and outer (£ = £ 2 ) boundaries as follows: 



(12) 



5+ i (a) = 0, n+(&)= 0, n-(6)=jo, (13) 
=0, ff- l (6) = 0, »_(&)= 0, (14) 

where i = 1, 2, • • •, m. We have totally 2m + 6 bound- 
ary conditions (0)-(|l^) for 2m+4 unknown functions 
$, En, n±, g±- ■ ■ g± m - Thus two extra boundary 
conditions must be compensated by making the posi- 
tions of the boundaries £1 and £ 2 be free. The two free 
boundaries appear because E\\ — is imposed at the 
both boundaries and because jo is externally imposed. 

2.4. TeV Fluxes 

The pair production cascade is primarily described 
by the curvature-radiated, GeV 7-rays. Therefore, 
the TeV emission due to Compton upscatterings is en- 
ergetically negligible and can be computed passively 
from the electrodynamic structure of the gap. Assum- 
ing that the azimuthal gap width is 0.5ttW/w-^c(= h) 
rad, we can evaluate the ratio between the TeV flux 
and the infrared one as (see §2.7 in Paper V) 
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where u>e refers to the solid angle in which photons 
with energy E are emitted; c^TcV = ^o.icv is assumed 
for simplicity. To estimate the upper bound, we take 
D± = TZ7LC- The inequality comes from the fact that 
the scattered photon energy cannot exceed the par- 
ticle energy. It is noteworthy that it is the infrared 
photons with energy ~ 0.1 eV that contribute most 
effectively as the target photons of comptonization. 
Neither the higher energy photons like surface black- 
body X-rays nor the lower energy photons like polar- 
cap radio emission contribute as the target photons, 
because they have either too small cross sections or 
too small energy transfer when they are scattered. 

Substituting the X-ray fields, £1 and /i of individual 
pulsars into the differential equations, we can solve 
and hence h; equation (15) then gives /r as a 
function of ai and jo. 

3. Application to Individual Pulsars 

In this section, we apply the theory to seven 
rotation-powered pulsars. We first present the ob- 
served X-ray properties of individual pulsars as input 
parameters in order of spin-down luminosity in ta- 
ble 1 . The X-ray spectrum of each pulsar is approxi- 
mated by separating into a soft blackbody component 
with temperature kT s and emitting area A s , a hard 
blackbody one with fcXh and , and a power-law one 
with dN x /de x = N p ie x a (e min < e x < e max ), where 
e min = 0.1keV/511keV and e max = 50keV/511keV 
are assumed for homogeneous discussion. 

Let us next consider the TeV fluxes. As a typical 
example, we present the results of /r for B1509-58 as 
a function of aj and jo i n figure [l|. The filled circles 
correspond to the maximum current density j cr above 
which no solution exists (see § 3.1 in Paper III and 
§ 3.1 in Hirotani & Okamoto 1998 for details) if all 
the other parameters are fixed. It follows from the 
figure that /r becomes maximum at jo — j cr (for a 
given ccj). We thus consider the case of jo = jcr in 
this section to estimate the upper bound of /r. 

The upper bound of /r for the seven pulsars are 
presented in table 2. For Crab and B0540-69, the syn- 
chrotron (or synchro-curvature) process dominates 
the simple curvature one when a; = 60°; therefore, 
this case is excluded in the table. For B1055-52, 
acceleration length scale for a particle to attain the 
saturated Lorentz factor exceeds the gap half width 
for ai — 60°. In another word, the mono-energetic 
approximation breaks down; therefore, the case of 



ct!j = 60° is excluded. It follows from the table that 
the TeV fluxes are unobservable (< 10 12 Jy Hz) for 
moderate infrared fluxes (< 10 9 Jy Hz), even when j 
is adjusted close to j cr . 



4. Discussion 

In summary, we have developed a one-dimensional 
model for an outer-gap accelerator. As the inclination 
angle (ai) increases, the gap center (where B z van- 
ishes) shifts starwards to be illuminated by a strong 
X-ray field, which reduces the pair production mean 
free path and hence the gap longitudinal width, W . 
At the same time, the starwardly shifted gap has a 
strong magnetic field, which increases the acceleration 
field (see eq. Q). As a result, the TeV flux increases 
with increasing a- u even though W decreases. How- 
ever, it is difficult to detect the TeV flux, even when a; 
is as large as 60°. We can therefore conclude that the 
difficulty of excessive TeV emission, which appears in 
the picture of Cheng, Ho, and Ruderman (1986a,b), 
does not arise in the present outer gap model. 

It is seemingly possible to obtain observable TeV 
fluxes, if ai exceeds 60°. In this case, the 'outer gap' 
is located close to the polar cap; the distance from 
the star is within 13% of otlc- The magnetic field 
there is so strong that the synchrotron (or synchro- 
curvature) process becomes important in general; this 
could be verified by computing the pitch angle evolu- 
tion of particles created in the gap (see §5.3 in Paper 
I). Such a 'polar-gap-like' outer gap will be discussed 
in a separate paper. 

For middle-age pulsars B0656+14 and Geminga, 
there is no stationary solution that satisfies the bound- 
ary conditions presented in § p.3j . However, if we al- 
low an external electric current flowing into the gap 
at the outer boundary, there exists a branch of so- 
lutions. In this case, the gap is no longer located 
around the null surface but shifts starwards. Interest- 
ingly, when the external current density approaches 
the Goldreich- Julian value, the 'outer gap' is found 
to shift to the polar cap and have a similar electro- 
dynamic structure of Scharlemann, Arons, & Fawley 
(1978). Such a unification of the outer gap and polar 
cap accelerators will be discussed in separate papers. 

This research owes much to the helpful comments 
of Dr. S. Shibata. The author wishes to express his 
gratitude to Drs. Y. Saito and A. Harding for valu- 
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Table 1: Input parameters and X-ray fields 



pulsar 


distance 


n 


logio M 


kT s 




kT h 


A h /A*l 


^V P i 


—a 


ref. 




kpc 


rad s _1 


lg(G cm 3 ) 


keV 




keV 




cm -3 






Crab 


2.49 


188.1 


30.53 










10 17.30 


1.8 


1 


B0540-69 


49.4 


124.7 


31.00 










10 14.15 


2.0 


2 


B1509-58 


4.40 


41.7 


31.19 










2Q14.04 


1.1 


2 


J1617-5055 


3.30 


90.6 


30.78 










10 12.64 


1.6 


3,4 


Vela 


0.50 


61.3 


30.53 


150 


0.066 










5 


B1951+32 


2.5 


159 


29.68 










2Q13.55 


1.6 


6 


B1055-52 


1.53 


31.9 


30.03 


68 


7.3 


320 


10 -3.64 
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t A* = AitrJ 2 represents the whole neutron star surface area; r* = 10 km is assumed. 

References: 1 Knight (1982) 2 Saito (1998) 3 Torii ct al. (1998) 4 Caswell et al. (1975) 

5 Ogelman et al. (1993) 6 Safi-Harb and Ogelman (1995) 7 Greiveldinger et al. (1996) 
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Fig. 1.— /r = {vF v ) TeV i '{vF v ) . leV vs. j for 
B1509-58, where {vF v )e indicates the vF v flux at 
photon energy E. The solid, dashed, and doted lines 
correspond to a\ = 30°, 45°, and 60°, respectively. 



Table 2: Expected 7-ray properties 



pulsar 


a ; 
deg 


jo = jcr 




/r 


t 


Crab 


30 


0.0485 


0.083 


6.1 x 


10 1 




45 


0.0371 


0.064 


1.8 x 


10 2 


B0540-69 


30 


0.0390 


0.074 


4.2 x 


10 1 




45 


0.0305 


0.053 


1.2 x 


10 2 


B1509-58 


30 


0.0375 


0.102 


5.1 x 


10° 




45 


0.0399 


0.083 


2.3 x 


10 1 




60 


0.0276 


0.052 


1.1 X 


10 2 


J1617-5055 


30 


0.0627 


0.136 


4.5 x 


10 1 




45 


0.0430 


0.104 


1.2 x 


10 2 




60 


0.0241 


0.049 


3.2 x 


10 2 


Vela 


30 


0.0620 


0.255 


2.8 x 


10 1 




45 


0.0375 


0.142 


4.2 x 


10 1 




60 


0.0264 


0.085 


1.7 x 


10 2 


B1951+32 


30 


0.0450 


0.133 


1.2 x 


10 1 




45 


0.0380 


0.095 


3.7 x 


10 1 




60 


0.0272 


0.057 


1.8 x 


10 2 


B1055-52 


30 


0.0343 


0.443 


2.1 x 


10° 




45 


0.0232 


0.182 


1.6 x 


10° 



t W indicates the full gap width. 

* h. = {fF v )Tev/(vF v )o.leV- 
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able advice. He also thanks the Astronomical Data 
Analysis Center of National Astronomical Observa- 
tory, Japan for the use of workstations. 
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